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Oceanic geostrophic turbulence is mostly forced at the surface, yet strong bottom-trapped 
flows are commonly observed along topographic anomalies. Here we consider the case of a 
freely evolving, initially surface-intensified velocity field above a topographic bump, and 
show that the self-organization into a bottom-trapped current can result from its tur- 
bulent dynamics. Using equilibrium statistical mechanics, we explain this phenomenon 
as the most probable outcome of turbulent stirring. We compute explicitly a class of 
solutions characterized by a linear relation between potential vorticity and streamfunc- 
tion, and predict when the bottom intensification is expected. Using direct numerical 
simulations, we provide an illustration of this phenomenon that agrees qualitatively with 
theory, although the ergodicity hypothesis is not strictly fulfilled. 



1. Introduction 

Bottom-intensified flows are commonly observed along topographic anomalies in the 
ocean. A striking example is given by the Zapiola anticyclone, a strong recirculation 
about 500 km wide that takes place above a sedimentary bump in the Argentine Basin, 
where bottom- i ntensified velocities of order 0.1 m.s ~'^ have been reported f rom in situ 
measurements ( Saunders Sz Kinslll995l ) and models (|de Miranda et al\\l99^ . This bot- 
tom intensifica tion may seem surprising , since the energy is mostly injected at the surface 
of the oceans ( Ferrari fc Wunschl[2009li . Indeed, the primary source of geostrophic tur- 
bulence is mostly baroclinic instability, extracting turbulent energy from the pot ential 
energy reservoir set at the basin scale by large scale wi nd patterns (iGiU et 
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involving surface-int ensified unstable modes, see e.g. (jSmithl l2007t iTulloch et al\ 12011 
IVenaille et al. 20116 ). It leads to the following question: how does a turbulent fiow driven 
by surface forcing evolve into a bottom-intensified fiow? A first necessary step to tackle 
this problem is to determine if an initially surface-intensified flow may evolve towards a 
bottom-intensified current without forcing and dissipation. Here we address this issue in 
the framework of freely-evolving stratified quasi-geostrophic turbulence, which allows for 
theoretical analysis with equilibrium statistical mechanics. 

The main idea of statistical mechanics is to p redict the large scale flow structure as 
the m ost probable outcome of turbulent stirring ( Miller et a/.lll992 : Robert fc Sommeria 



1992| ). The interest of the approach lies in the possibility to compute this large scale 



flow structure and to study its p roperties with respect to a few key parameters given 
by dynamical invariants, see e.g. (|Maida fc Wan3l2006l ; iBouchet fc VenaiUe 2011.) . We 
focus here on a particular class of equilibrium states amenable to analytical treatment, 
namely, those characterized by a linear relation between potential vorticity (PV) and 
streamfunction. 
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Phenomenol oKical argume nts for the formation of bottom-trapped flows were previ 
ously given by IPewai (1998) in a forced-dissipat ed case. We provi de here a complemen 
tary point of view, built upon previous results of Merrvfield 1 19981 ) . who computed criti 
cal states of equilibrium statistical mechanics for truncated dynamics. He observed that 
some o f these states were bottom intensified in the presence of topographv. I Venaille et "al 



I 2011cl ) showed how to find the equilibrium states among these critical states, and how 
they depend on the initial fine-grained enstrophy profile. However, they focused on the 
role of the beta effect in barotropization processes and did not discuss the effect of bottom 
topography. Here we show that bottom-trapped currents are actual statistical equilibria 
in the presence of sufficiently large bottom topography, and we present direct numeri- 
cal simulations of the free evolution of an initial surface-intensified flow towards these 
bottom-trapped currents. 



2. Equilibrium states of quasi-geostrophic flows 

Here we consider an initial value problem for a Boussinesq, continuously stratified 
quasi-geostrophic fluid. Such flows are stably stratified with a prescribed buoyancy profile 
N(z) above a topography anomaly hb(x, y), and are strongly rotating at a rate /o/2 such 
that UL/ fo ^ 1, where U and L respectively are typical horizontal velocity and length of 
the flow. At leading order, the flow is at geostrophic balance: horizontal pressure gradients 
compensat e the C oriolis term. The dynamics is then given by quasi-geostrophic equations 
(see Yalli3 (|2006h section 5.4): 



dtq + d^ipdyq - dyipd^q = 



(2.1) 



with dzi^lz 



d_ 

dz 

= 0, 



111 

7V2 dz 



A 



-hh 



(2.2) 



(2.3) 



z=-H 



where ip is the streamfunction, A the horizontal Laplacian, H the ocean depth. We 
consider Cartesian coordinates, but the earth sphericity is taken into account at lowest 
order through the term /3j/, which is due to the meridional variations of the projection of 
the earth rotation rate on the local vertical axis. This is the beta plane approximation. 
The lateral buoyancy variations have been neglected in the upper and lower boundary 
conditions. Quasi-geostrophic approximation neglects topography variations to leading 
order. This is why the term hi, appears only at the boundary z — ~H . 

The dynamics (| 2. 1112.2112731) takes place in three dimensions, but is expressed as the 
advection of a scalar field (the PV q) by a non-divergent horizontal velocity field, as 
in two-dimensional Euler dynamics. For this reason, the phenomenology of stratified 
quasi-geostrophic turbulence is similar to two-dimensional turbulence: there is an inverse 
energy cascade, self-organization at the domain scale, and no dissipative anomaly. 

Quasi-geostrophic dynamics develops complex PV filament s at finer and finer scales at 
each depth z, associated with a forward cascade of enstrophy (jKraichnan fc Montgomery 
198d) . Rather than describing the fine-grained structures, equilibrium statistical theories 
of two-dimensional turbulent flows, by assuming ergodicity (or at least sufficient mixing 
in phase space, corresponding to stir ring in physical space), predict self-organizatio n 
of the flow on a coarse-grained level ( Robert fc Sommerial (Il992l ): iMiller et al\ (|l992l ). 
MRS hereafter): a "mixing entropy" is maximized by taking into account the constraints. 
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namely, all the flow invariants, which are the total energy 



dz / dxdy 
H Jv 



(2.4) 



and the Casimir functionals Cg{z)[q\ — dxdy g{q) where g is any continuous function, 
and T) the domain where the flow takes place. Here we consider the case of a doubly- 
periodic domain. In the case g{q) = q'^ /2, Cg is the enstrophy. The output of the theory 
is a coarse-grained PV field obtained by locally averaging the fine-grained PV field q. 

Let us call Eq = £{qo) and Zo{z) — (1/2) J^dxdy q^ the energy and fine-grained 
en strophy, r espectively, of the initi al condition given by (? o(x, y, z). Using a general result 
of Bouchet (|2008i) . it is shown in Venaille et all (|2011cl) . Appendix 1, that in the low 



energy limits the calculation of MRS equilibrium states amounts to finding the minimizer 
Qmin of the "total coarse-grained enstrophy" 

1 



dz 



-H 



dxdy 



among all the fields q satisfying the energy constraint 



= \ i da;dy fohti'lz^-H - \ ( dxdy / dz {q 

^ Jv ^ Jv J ~H 



py)^^Eo, 



(2.5) 



(2.6) 



where Eq. (|2.6p is obtained after integrating Eq. (|2.4p by parts. Independent of the 
statistical mechanics arguments, the variational problem (|2.5ll2.6p can be seen as a gen- 
er alization to the stratified case o f the phenomenological minimum enstrophy principle 

of Bretherton fc Haidvoeell (1197a). Besides, these states are statistical equilibria of the 

^ " — — ' ' II — — k 

truncated dynamics (jKraichnan fc Montgom erv 1980; Sal mon et al\\197w in the limit of 

infinite wave- number cut-off ( Carnevale fc Freder iksen 198?!) . 

Critical states of the variational problem (|2.5ll2.6p are computed by introducing La- 
grange multiplieJil /3t associated with the energy constraint, and by solving SZ^p""^ + 
j3t5£ = 0, where variations of the functionals are taken with respect to q^ leading to the 
linear relation q = /SfZoip, see also I Venaille et all (j2011cl ). The next step is to find which 
of these critical states are actual minimizers of the coarse-grained enstrophy for a given 
energy. 



3. Analytical calculations and simulations 

We consider the configuration of Fig. [T]-a: the stratification is linear in the bulk {N is 
constant for —H + h < z < —h), and homogeneous in two layers of thickness h <^ H at 
the top and at the bottom, where iV = 0+. In these upper and lower layers, the stream- 
function is depth independent, denoted by iJ;'^°p {x , y , t) and {x , y , t) , respectively. In 
these layers, the dynamics is then fully described by the advection of the vertical av- 
erage of the PV fields, denoted by q*°P{x,y,t) and q''"*' {x , y , t) . The interior PV field is 



denoted by g™*(x, y, z 
by inverting the following equations 



t). For a given field q*°P, q™*, g''°*. 



JLI. 

hN^ dz 



the streamfunction is obtained 



(3.1) 



^bot 



/3y-/o^=AV. 



hot 



hN^ dz 



(3.2) 



=h-H 



f The inverse temperature /3t should not be confused with the /3 eflfect. 
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Figure 1 . a) Sketch of the flow configuration. The continuous red hne represents the initial fine- 
grained enstrophy profile (here Z™* = 0, as assumed in the Appendix). The dashed blue line 
represents the density profile, and the dashed-dotted green line represents the streamfunction 
amplitude shape, which is initially surface-intensified. The thick continuous black line repre- 
sents bottom topography, b) Vertical slice of the meridional velocity field v of the statistical 
equilibrium state in the low energy limit. 



g"* -/3y = 



-ip ior - H + h < z < -h 



(3.3) 



11^ 

= i:{x, y, -h), ^""^ = i:{x, y, ^H + h). (3.4) 

Equations (|3.HI3.2p are obtained by averaging Eq. (|2.2[) in the vertical direction in the 
upper and the lower layers, respectively, and by using the boundary condition (|2.3p . 
In the following, the initial condition is a surface-intensified velocity field induced by a 
perturbation of the PV field confined in the upper layer: 



pert 

% 



bot ^ foj 



(3.5) 



It is assumed in the following that /3y <C q^"^^* ^ fohb/h. The PV fields are therefore as- 
sociated with fine-grained enstrophies Zq"* <^ Zq"^ <^ Zq°* . We compute in the Appendix 
the coarse-grained enstrophy minimizers of this configuration by solving the variational 
problem (|2.5ll2.6p . The main result is that for a fixed topography, in the low energy limits 
the equilibrium streamfunction is a bottom-intensified quasi-geostrophic flow such that 
bottom streamlines follow contours of topography with positive correlations, see Fig. [1]- 
b. Given the choice of our initial condition, considering the low energy limit for a fixed 
topography is equivalent to considering a large topography limit for a fixed energy. 

We perform simulations of the dynamics by considering a horizontal discretization of 
512^ in a doubly periodic square domain of size L — 27r, and a vertical discretization 
with 10 layers of equ al depth h = 0.1, using a pseudo-spectral quasi-geostrophic model 
( Smith fc Vallij 2001 ). The linearly stratified region of Fig. [T]is therefore approximated 
by 8 layers of depth h. We present results without small scale dissipation, but we checked 
that adding a weak enstrophy filter or hyperviscosity did not change qualitatively the 
results. 

The perturbation Qq'^'^* of the PV field in the upper layer has random phases in spectral 
space with a Gaussian power spectrum peaked at wavenumber Kq = 5 , with variance 
6Ko = 2, and normalized such that the total energy is equal to one {Eq = 1), the latter 
amounting to a rescaling of time. The bottom topography hb is a Gaussian bump of 
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Figure 2. a),c),e) PV field at three successive times. In each panel, only layers 1 (top), 5 
(middle) and 10 (bottom) are represented. b),d),f) Vertical slices of the meridional velocity 
fields V taken at the center of the domain (j/ = 0), and associated with the PV fields given in 
panels a),c),e), respectively. The bold continuous dark line represents bottom topography. The 
continuous black contours of panel f) give the structure of the statistical equilibrium state in 
the low energy (or large topography) limit, corresponding to Fig. 1-b. Contour intervals are the 
same as those between the different shades. 

amplitude Hh and of typical width a = L/2. There are three independent parameters: 

1 /2 

where Vq — {Eq/L^H^ is an estimate for the velocity. The first parameter is the ratio 
between the first baroclinic radius of deformation and the domain scale, taken to be 
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Ld/L — 1/10. The second one is the ratio between the Rhines scale and the domain 
scale. We choose Lp/L — 1/3. The third one is the ratio of a topographic Rhines scale 
to the domain scale, taken to be Lh^/L = 1/10. We show in the Appendix that the 
latter choice is consistent with the low energy or large topography limit considered in 
the previous section. The choices of /3 and Hi, are also consistent with the hypothesis 
that Zg"* <^ Zg°^ <C Zq°^. The choice Kq = 6 for the initial condition is motivated 
by the fact that in the forced-dissipated case, the energy is injected through baroclinic 
instability with a most unstable mode having a w avelength of a few Rossby wavelength 
of deformation, see e.g. ()Smithl 120071: IVallidl200d ). We checked that changing this mitial 
horizontal length scale did not change the results. 

We s ee in Fig.[2l - b that the initial condition resembles a classical surface quasi-geostrophic 
mode ( Held et al\\l993i) . After a few eddy turnover times, the enstrophy of the upper 
layer has cascaded towards small scales as shown by numerous filaments in Fig. [2]-c, 
concomitantly with an increase of the horizontal energy length scale. The projection of 
a surface quasi-geostrophic mode on a Fourier mode of wavenumber K is characterized 
by a typical e-folding depth of fo/NK on the vertical. The inverse energy cascade on 
the horizontal leads therefore to a deeper penetration of the velocity field, shown in Fig. 
[2]-d. When this velocity field reaches the bottom layer, it starts to stir the bottom PV 
field. This induces a bottom-intensified flow, which then stirs the surface PV field, and 
so on. The flow reaches a stationary state after a few dozen of eddy turnover times: the 
q — tp relation in each layer converges towards a well defined functional relation after 
coarse-graining. The existence of this functional relation implies dtq — 0, meaning that 
there are no temporal fluctuations. The corresponding flow is shown in Fig. [2]-e, which 
clearly represents a bottom-intensified anticyclonic flow above the topographic anomaly, 
qualitatively similar to the one predicted by statistical mechanics in the low energy limit. 
How ever, we observe the form ation of multiple jets (here 2) around the topographic bump, 
as m IVallis & Maltrildl which indicates that the "sufficient stirring" hypothesis is 

not strictly fulfilled in the simulations. 

The same evolution from a surface-intensified turbulent field towards a bottom-trapped 
flow was observed in the absence of a /3-effect, except that there remained two small 
surface-intensified coherent vortices locked above the topograph y extrema. T hese vortices 
are destroyed by radiation of Rossby waves in presence of /?, as in iMaltrud fc Vallis ( 199lh . 
In a way, small values of /3 enhance PV stirring, even if they have a negligible effect on 
the structure of the equilibrium state. However, if /? is too large, then Lp <^ Lh^, which 
leads to zonal jets that prevent the formation of bottom-trapped flow along topographic 
anomalies. 



4. Conclusion 

For sufficiently large values of bottom topography, statistical mechanics predicts that 
an initially surface-intensified velocity field will evolve towards a bottom-intensified fiow 
following contours of topography. To the best of our knowledge, this is the first analyti- 
cal calculation of statistical equilibrium states for stratified quasi-geostrophic equations 
abo ve topo graphy, although there have been many works on one or few layer mod e ls, see 
e.g. Maida fc Wand |2006f ). and also, critical states were computed bv lMerrvfieldl ( 19981 ) 



in the stratified case. We chose to present these calculations in the simplest possible set- 
ting in order to give the main physical ideas. We found qualitative agreement between 
statistical mechanics predictions and a direct numerical simulation of the freely evolving 
dynamics. 

The statistical mechanics approach highlights the crucial role of energy and fine-grained 
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enstrophy conservation (among other invariants) while taking into account the turbulent 
nature of the dynamics. Here we considered the equilibrium states characterized by energy 
and enstrophy only, which is justified in a low energy limit, and which leads to linear q — ijj 
relations. The next step will be to consider higher order invariants, i n order to account for 
equilib rium states associated with non-linear q—ijj relations, see e.g. Bouchet fc Simonnetj 
in the context of Euler equations. 



The main shortcoming of the statistical mechanics approach is that it predicts only 
the final state organization, and not how the system evolves towards this state. In ad- 
dition, there remain quantitative discrepancies between theory and simulations, because 
the ergodicity, or "sufRcient stirring" hypothesis may not be strictly fulfilled in some 
cases. Addressing specifically the range of parameters for which the statistical mechanics 
predictions are accurate, and estimating how the relaxation time towards equilibrium 
depends on these parameters will be the object of a future work. 

These results have important consequences for ocean energetics: topographic anoma- 
lies allow transferring surface-intensified eddy kinetic energy into bottom-trapped mean 
kinetic energy, which would eventually be dissipated in th e presence of bottom friction , 



as for instance in the case of the Zapiola anticyclone (jDewar.,1998 : Venaille et al\\20llc$ ). 
In the case of the Zapiola anticyclone, the dissipation time scale is of the order of a few 
eddy turnover-time. It is therefore not a priori obvious that the results obtained in a 
freely evolving configuration may apply to this case. However, we consider that it was a 
necessary first step. One can now build upon these results to address the role of forcing 
and dissipation in vertical energy transfers above topographic anomalies. 
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5. Appendix: calculation of the solution 

We compute statistical equilibria for the configuration sket ched in F ig. [J and detailed 
in section 3, by applying a general method presented in .Venaille fc Bouc het (2011,). Here 
we consider the limit of smal l beta effect. The cas e without topography but including 
the beta effect is discussed in I Venaille et all ( 2011cf ). 

When /3 — ^ 0, only the upper and lower layers are dynamically active. The reason is 
that the interior fine-grained enstrophy tends to zero. We have seen in section 2 that 
a solution of the variational problem (|2.5ll2.6p necessary satisfies q = PtZoi/j. This gives 
Jjjdxdy q'^/Zo ^ Zq. We conclude that lim2o_i.o+ /p da;d?/ if/Zo — 0. The variational 
problem ()2.5ll2.6p amounts therefore to finding the minimizers of 

ld.dy + ld.dy (q'-'f, (5.1) 

with the constraint 

One could also derive this variational problem directly from the Miller- Robert-Sommeria 
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variational problem, using the method presented in the Appendix o f IVenaille efd\ \20U(h . 



and assuming that PV levels are zero everywhere except in the the lower and upper ho- 
mogeneous layers. 

In order to solve the variational problem (|5.m5.2p , we use the fact that any minimizer 
^/iT' ^A* ^^^"^ "free energy" T — Zl°* + /3t£ is a minimizer of the enstrophy with the 
energy constraint Ep^ = £ , see e.g. ( Venaille fc Bouche^l2011 ) and references 

therein. 

The first step is to compute the free energy. For that purpose, it is convenient to 
consider Fourier modes ij) = '^f, ^ipk,iexp {ikx + ily) and solve Eq. p.3ll3.4p for each 
Fourier component: 

,A Ttop sinh^£-h^7A^ 2bot sinh(z/Ajf) fa . 2 , ,2n1/2 

3i,h(iJ/A^) -^'^■'sinhWA;,) ' ^^' = ]^'^=(^ +0 ' 

(5.3) 

where we considered the limit h H to simplify the expressions. Each Fourier component 
of the streamfunction is a linear combination of a surface quasi-geostrophic mode and 
a bottom-trapped quasi-geostrophic mode. Substitution of Eq. (j5.3p into the Fourier 
projections of Eqs. (|3. 1113. 21) yields 

with = -X^ ^ + — -T^TTT-T V b,=K---^^^. (5.5) 



htim\i{H/\K) J ' hsmh{H/\K 
In order to compute the free energy F = -I- /3tf , we first express the enstrophy (|5.ip 
and the energy ()5.2p in term of the Fourier components of the streamfunction, potential 
vorticity fields, bottom topography. We finally use Eq. (|5.4p to express the streamfunction 
in term of the potential vorticity fields and the topography, which yields 





2/3* c*op ( c^ot 

2 IT^k.l I 1k,l ' 

After some manipulations, we find that the quadratic part of J- is positive definite for 
/3t e (/3™™,+oo), where /Sj"™ is the largest root of (/^tZ*"^ - ai) (/3t^^°* - ai) = h\, 
implying that there exits a unique minimizer for each /3t £ (/3™™, -l-oo) . In addition, 
there is at least one unstable (negative) direction for T when /3t G (—00, /?™*") , implying 
that there is no minimizer for this range of parameters. 

When it exists, the minimizer g„jj„ of the functional J- can be explicitly computed 
for each value of /3f. The first step is to compute the critical points of the functional, 
which are solutions of 8T = SZ^g* + PtSE = 0, where variations are taken with respect to 
the PV fields fP. This leads to = /^t^oV'* with i e {bot,top}. Substituting the 
Fourier projections of these expressions in Eq. ()5.4p yields 

^top _ fo bK /r \ g-bot _ PtZl°^ ~ aK -rtop 

- h ^p,zT - aK) {f^tZr - aK) - h\ M ' - bK ^'^'^ ■ 

(5.7) 

For a given non-zero topography, the energy tends to zero when /3t — -t-00, and tends 
to infinity when /3t — >■ /?""". It is therefore possible to reach any admissible energy by 
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varying the inverse temperature between /3™™ and +00. It thus means that we have 
found all the energy-enstrophy statistical equilibria. 

Let us first consider the limit /3fZQ°*i^ — > +00. The solution (|5.7p tends to 



fobx (hb) fo (hh) 

\ / k,l ^hnt ^ ' k,l 



We see that ■0''°* = fohf,/ {h/3tZQ°*) is along contours of topography, and that the vertical 
structure of the equilibrium state is dominated by the bottom-trapped quasi-geostrophic 
mode, since <€. ijj^"*- Using the energy expression (|2.4p . we obtain Ea ^ HL"^ (^foHb/Lh(3tZQ° 
and substituting this expression in Eq. (|32^iii) yields /StZoL'^ - {L/LhJ . These esti- 
mates explain a posteriori the term low energy limit or large topography limit when 
/3tZQ°*L^ — -l-cxo, and the fact that this limit corresponds to Lh,^ ^ L. 

In the low topography limit or large energy limit, the ratio Zq°*/Zq°^ tends to zero, 
and we deduce from Eq. (|5.7p that f3t ~ /3™™ « (ai — 6f /ai) /Zq°^, which implies that 
^top jg dominated by its Fourier component on the horizontal mode K = 1. We also 
deduce from Eq. (|5.7p that w — 6i/aiV'*°^ ^ -0*°^: the equilibrium state is a surface 
quasi-geostrophic mode. 



We have thus described a class of equilibrium states varying from a surface quasi- 
geostrophic mode condensed in the gravest Laplacian eigenmode {low topography or large 
energy limit) to a bottom-trapped quasi-geostrophic mode following contours of topog- 
raphy {low energy or large topography limit). In the framework of MRS theory, only the 
low energy limit is consistent with the derivation of the quadratic variational problem 
(lOIOl . 

Finally, the numerical simulations were performed with ten homogeneous layers, which 
is a slightly different case from the one studied in this Appendix, with linear interior 
stratification. However, the present analytical results easily generalize to arbitrary sta- 
ble interior stratification profiles. We do not expect different qualitative behaviors for 
the equilibrium states, although the convergence towards equilibrium may depend on 
stratification properties. 
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